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ABSTRACT 



00 ! iZST/NICMOS PSF-subtracted coronagraphic observations of HD 181327 

have revealed the presence of a ring-like disk of circumstellar debris seen in 1.1 
(Urn light scattered by the disk grains, surrounded by a diffuse outer region of 
lower surface brightness. The annular disk appears to be inclined by 31° 7 ± 1?6 
from face on with the disk major axis PA at 107° ± 2°. The total 1.1 /im flux 
O ' density of the light scattered by the disk (at l'/2 < r < 5'/0) of 9.6 mJy ± 0.8 mJy 



is 0.17% ± 0.015% of the starlight. Seventy percent of the light from the scatter- 
ing grains appears to be confined in a 36 AU wide annulus centered on the peak 
of the radial surface brightness (SB) profile 86.3 AU ± 3.9 AU from the star, well 
beyond the characteristic radius of thermal emission estimated from IRAS and 
Spitzer flux densities assuming blackbody grains (~ 22 AU). The 1.1 /im light 
scattered by the ring: (a) appears bilaterally symmetric, (b) exhibits direction- 
ally preferential scattering well represented by a Henyey-Greenstein scattering 
^ ■ phase function with g HG = 0.30 ± 0.03, and (c) has a median SB (over all az- 

imuth angles) at the 86.3 AU radius of peak SB of 1.00 ± 0.07 mJy arcsec -2 . No 
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photocentric offset is seen in the ring relative to the position of the central star. 
A low surface brightness diffuse halo is seen in the NICMOS image to a distance 
of ~ 4". Deeper 0.6 /zm HST/ ACS PSF-subtracted coronagraphic observations 
reveal a faint (V ~ 21.5 mag arcsec~ 2 ) outer nebulosity from 4" < r < 9", asym- 
metrically brighter to the North of the star. We discuss models of the disk and 
properties of its grains, from which we infer a maximum vertical scale height of 
4-8 AU at the 87.6 AU radius of maximum surface density, and a total maximum 
dust mass of collisionally replenished grains with minimum grain sizes of ~ 1 /im 

Of « 4 M M oon- 

Subject headings: circumstellar matter — infrared: stars — planetary systems: 
proto-planetary disks — stars : individual (HD 181327) 

1. INTRODUCTION 

Fifteen years elapsed between the acquisition of the prototypical scattered-light image 
of a debris disk around a young main-sequence star ([3 Pictoris; Smith & Terrille 1984), 
offering strong support to the conjecture of the "Vega phenomenon" (far-IR emission in 
excess of a stellar photosphere from orbiting thermally emissive dust (Backman & Paresce 
1993) as suggested from earlier IRAS observations (Aumann et al. 1984)), and the second 
example of a dusty circumstellar disk - in the form of a ring - seen by starlight illuminating 
its constituent grains (HR 4796A; Schneider et al. 1999). Advances in space- and ground- 
based high contrast imaging capabilities have given rise to increasing, though still very small, 
numbers of scattered-light images of circumstellar debris systems. The very recent discoveries 
of a high surface brightness (3 Pictoris-like nearly edge-on disk around the A-star HD 32297 
(Schneider, Silverstone & Hines 2005 [=SSH05]), and the very low surface brightness 140 AU 
ring around an A-star much older than f3 Pictoris, i.e., Fomalhaut (Kalas, et al. 2005) speak 
to the diversity in disk morphologies and scattering properties. While ring-like structures are 
obscured in scattered-light images of nearly edge-on disks, their intrinsic morphologies may 
be inferred by inversion of their radial surface brightness profiles (egs.: Buitrago & Mediavilla 
(1986), Pantin & Lagage (1996), Augereau et al. (2001) for j3 Pictoris; and models presented 
by Krist et al. (2005) and Metchev et al. (2005) for AU Microscopii) . We now report the 
discovery of a ring-like disk of circumstellar debris seen about the star HD 181327. This new 
addition to the growing menagerie of circumstellar debris disks is, in many respects, similar 
to the HR 4796A star:disk system. The recent imagery of a broader ring about the older 
solar analog, HD 107146 (Ardila et al. 2004), suggests that such structures are a common 
(if not prevalent) outcome of the planet formation process. 
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HD 181327 (HIP 95270) is an F5/F6V star (T eff « 6450 K; Nordstrom et al. (2004), 
log(g) = 4.5) at 50.6 pc ± 2.1 pc (Perryman et al. 1997) with a Vega-like far-IR excess noted 
by Backman & Paresce (1993). HD 181327 was identified as a main sequence debris-disk 
candidate by Manning & Barlow (1998) with an L ir / L* = 0.2%. Kalas et al (2002) favored 
the Vega phenomenon interpretation of its far-IR excess based upon its youth and its gas 
depleted environment. The star's youth was first suggested from its potential membership 
in the ~ 40 Myr Tucanae association (Zuckerman & Webb 2000), but later identified as 
a member of the younger (~ 12 Myr) (3 Pictoris moving group (Zuckerman & Song 2004; 
Mamajek et al. 2004). 

The conjecture of a stellar age of ~ 10 7 yr is further supported by two independent 
observations. First, HD 181327 appears to exhibit common proper motion with the 7' (0.1 
pc at 50 pc) distant HR 7329, also identified as a member of the (3 Pictoris moving group 
by Zuckerman et al. (2001). Both stars have statistically identical Hipparcos parallaxes 
and (long baseline) Tycho 2 proper motions. HR 7329 is a ZAMS A-star possessing a close 
substellar companion 1 with an isochronal age of ~ 10-30 Myr (Lowrance et al. 2000). Second, 
HD 181327's Hipparcos derived B-V color index of +0.480 is completely consistent with an 
unreddeded F5-F6 star 2 . This, then, is in agreement with HD181327's e(b-y) = 0.002 mag 
(Nordstrom et al. 2004), and the stellar luminosity is then 3.1 ± 0.3 L sun from which the 
stellar age estimated using evolutionary tracks by D'Antona & Mazzitelli (1997), Palla & 
Stahler (2001), and Baraffe et al. (1998) is consistent with ~ 20 Myr. 

Assuming blackbody grains, Zuckerman & Song (2004) suggested a thermal equilibrium 
temperature of 65 K for the IR emissive dust at a characteristic orbital radius 35 AU from 
HD 181327. Neuhauser et al. (2003), using NNT + SHARP-I & SofI, and Chauvin et al. 
(2003), with ESO 3.6 m AO + ADONIS/SHARP-II, conducted searches for a companion 
candidate to HD 181327, but found none. 



lr The total mass of HD181327 + HR7329A/B is ~ 3.8 M Sm , and stellar binaries with separations ~ 0.1 
pc are not unheard of with similar systemic masses. Hence, in light of recent discussion on brown dwarf 
formation and binary populations (egs., see Close et al. (2003) fig 15 and Burgasser et al. (2006)) one might 
conjecturally consider HD181327 + HR7329A/B to be a triple system. 

2 In the Hipparcos catalog, the F5 and F6 stars within 75 pc with parallaxes of S/N > 8 have median 
B-V color indices of +0.464 and+0.496, respectively. Hence, there appears to be no evidence for reddening 
for HD 181327 (unsurprising given its distance and environment). 
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2. NICMOS OBSERVATIONS & DATA PROCESSING 

During HST cycle 13 we conducted a NICMOS coronagraphic imaging survey (HST 
GO/10177; PI: Schneider) of 26, ~ 10 Myr and older, main sequence stars with strong far-IR 
excess emission, including HD 181327, to identify and provide spatially resolved scattered- 
light images of previously suspected, but unimaged, dusty circumstellar disks (SSH05). 
HD 181327 was observed at two field orientations in a single spacecraft orbit on 2005 May 02 
UT using HST/NICMOS camera 2 (scale: ~ 75.8 mas pixel" 1 ) with F110W (A e// = 1.104 
/im, FWHM = 0.5915 /im) filter direct and coronagraphic imaging, and supporting F180M 
(X e ff — 1-7968 /zm, FWHM = 0.0684 //m) target acquisition exposures. The specifics of the 
observations, and those of the PSF reference stars from which our PSF-subtracted imagery 
we discuss in §2.1, are given in Table 1. Details of the observing strategy, data reduction 
and calibration methodologies, and image combination processes which we apply to all de- 
bris disk candidate (and reference PSF) star observations in our HST GO/10177 survey are 
discussed by SSH05 and to which we refer the reader for details. 



2.1. NICMOS PSF-Subtractions 

To further reduce the level of the residual instrumentally diffracted and scattered cir- 
cumstellar light, beyond the improvement in image contrast afforded by the coronagraph 
itself, we subtract astrometrically registered, flux scaled images of coronagraphic reference 
PSFs of non-disk stars of similar spectral energy distributions from each of the two (one 
at each field orientation) calibrated combined images of HD 181327. We follow the pre- 
cepts of the two- orient at ion coronagraphic PSF-subtraction process described by SSH05. 
We selected ten candidate PSF stars with near-IR colors similar to HD 181327. Two of 
these were specifically programmed very bright PSF calibration stars with no IRAS detected 
excess emission, HR 4748 (J=5.52 mag) and HIP 7345 (J=5.49 mag), each observed at a 
single field orientation using a full spacecraft orbit and yielding very high SNR. Eight other 
PSF template candidates were selected from GO 10177 survey targets that exhibited no 
evidence of scattered-light excesses from circumstellar dust. Each of these eight stars were 
observed at two field orientations (thus providing two PSF templates each), but were much 
less deeply exposed than our PSF calibration targets due to their lesser brightness and the 
detailed nature of the observing/integration time strategy for these targets. All 36 HD 
181327-minus-PSF template subtractions revealed an elliptical ring-like disk that rotated in 
the celestial frame, following the 29?9 difference in target orientation between the two HD 
181327 observational visits (e.g., Fig 1A-D). PSF-minus-PSF subtractions exhibited no such 
structures (e.g., Fig IE) while in all cases rotationally invariant features due to instrumental 
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scattering and diffraction remained fixed in the frame of the detector. 

The four images resulting from subtracting the HR 4748 and HIP 4375 high SNR refer- 
ence PSFs from each of the two images of HD 181327 (differing in orientation 29?9) were of 
higher photometric efficacy than those obtained using the serendipitous "ensemble" of PSFs 
derived from the GO 10177 survey targets and, generally were more devoid of azimuthally 
anisotropic high and mid spatial frequency artifacts arising from aliasing in imperfect PSF- 
subtractions. The HD 181327 debris disk image, which we discuss in this paper, was created 
using these four PSF-subtracted images. The four images were median combined with equal 
weight after alignment, rotation to a common orientation, and masking regions in individual 
images affected by diffraction spikes from the HST secondary mirror support structure (see 
Figure IF). The final NICMOS image, shown in Figure IF (and 2A), and discussed in §§ 5 
& 6, has a total integration time of 1536 s in most pixels. 

3. ACS OBSERVATIONS 

HD 181327 was observed earlier in HST/GTO program 9987 (PI: Ford) on 2003 Novem- 
ber 02 UT using the HST Advanced Camera for Surveys/High Resolution Camera (ACS/HRC; 
scale: ~ 25 mas pixel -1 ) in direct and coronagraphic imaging modes. Two 5 s exposures 
were taken without the coronagraph in the F606W ("wide" V-band) filter. After a 0.5 s 
F502N acquisition exposure, the star was placed behind the 0'.'9 radius coronagraphic occult- 
ing spot and two 1235 s F606W exposures were taken. The ACS images were calibrated by 
the HST OPUS 3 pipeline, and the two coronagraphic exposures were combined with cosmic 
ray rejection. 

3.1. ACS Coronagraphy & PSF Subtraction 

By suppressing the diffraction pattern of the occulted star, the ACS coronagraph reduces 
the surface brightness of the wings of the stellar PSF by about a factor of 7 - 10. As 
with NICMOS, an additional gain in image contrast can be achieved by subtracting the 
underlying stellar PSF using a coronagraphic image of another occulted star of similar color 
reducing the remaining "background" by factors of 10 - 100. Additional information on ACS 
coronagraphic performance is presented in § 6.2 of Gonzaga et al. (2005). 

In the case of HD 181327, observations of the F5V star HD 191089 were used to create 



3 lit t p : / / www .stsci.edu/ resources /softwareJiardware / opus / 
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a template PSF image for subtraction. HD 191089 was imaged in a single spacecraft orbit 
immediately following the observations of HD 181327, and a very similar imaging and data 
processing strategy was employed. Two 6 s direct and two 1005 s coronagraphic (r = 0'. / 9 
occulting spot) exposures were obtained with the F606W filter. The HD 191089 PSF tem- 
plate was scaled to the flux density of HD 181327 4 . After finite-aperture photometric and 
geometrical distortion corrections, the target-to-PSF scaling obtained agreed within 2% of 
expected based upon the V-band magnitudes of the stars. The subtraction of HD 181327 
was optimized by iteratively shifting and subtracting the flux-scaled HD 191089 image until 
residuals induced by mis-registration, beyond 2" from the occulted star, were minimized and 
spatially symmetric. This process co-registers the stars to ~ 0.05 HRC pixels (1.25 mas). 
The final subtracted image (shown in figure 3 at three different image scales and display 
stretches) was corrected for the well-characterized geometric distortion in the HRC camera. 



3.2. ACS Image Data & Caveats 

Near the coronagraphic obscuration, a residual halo of starlight unsuppressed by the 
ACS coronagraph arises from instrumental scattering (seen as radial streaks) and diffraction 
by the occulter itself (seen as concentric rings centered on the occulter). The scattered-light 
streaks are relatively insensitive to star-to-occulter positioning, but the ring-like diffracted 
light distribution is sensitive to the location of the star within the footprint of the occul- 
ter. The ACS PSF-subtracted image of HD 181327-minus-HD 191089 (after flux-scaling 
and image registration) shows instrumentally-induced residuals beyond the edge of the oc- 
culting spot. The residuals are dominated by diffraction rings arising, predominantly, from 
differences in the positioning of the two stars behind the occulter. The registration of the 
HD 181327 and HD 191089 images required a larger than typical translation (30.1 mas) of 
the reference PSF indicative of a target centration error in the acquisition process (assuming 
the position of the occulting spot remained stable over the short period of the two-orbit 
sequential observations of both stars). The position error was much greater along the image 
horizontal (w.r.t. the detector), so the PSF-subtraction residuals are also larger along that 
direction. Due to the asymmetry of the PSF subtraction- residuals near the edge of the spot: 
(a) the HD 181327 disk can be detected reliably from I'M < r < 2"0 from the star only 



4 To establish the flux density scaling, for both stars the number of saturated pixels in the PSF core 
was multiplied by 165000 e~ pixel -1 (full well depth of an HRC CCD pixel) and added to the flux density 
measured within a 1" radius aperture. Charge in excess of a pixel full-well is transferred, without significant 
inefficiency, to neighboring non-saturated pixels. For additional information see §§ 4.2.6 and 4.3.1 of Gonzaga 
et al. (2005), Gilliland (2004) and Bohlin et al. (2003). 
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within ~ ± 10° arcs centered on the image-Y axis on opposite sides of the star (figure 3A); 
(b) from 2" < r < 4" a diffuse halo is seen (figure 3B), though its near circular-symmetry and 
ring-like substructures might implicate a residual contribution of instrumental origin (from 
imperfect PSF-subtraction) that is somewhat difficult to quantify. 

Further from the coronagraphic occulting spot, from 4" < r < 9" (well beyond the 
distance leading to PSF-subtraction artifacts discussed in (a) and (b)), within ± 45° of 
North a faint nebulosity with a mean surface brightness of V ~ 21.5 mag arcsec -2 is seen 
(figure 3C). 

4. Spitzer/MIPS and IRAS PHOTOMETRY 

New 24, 70, and 160 /xm imaging photometric observations of HD 181327 were obtained 
with the MIPS instrument on the Spitzer Space Telescope on 2004 April 06 in GTO program 
72 (F. Low, PI). We use the photometric measures derived from these observations, along 
with IRAS Faint Source Catalog (FSC) 25, 60, and 100 fim photometry, to better define the 
thermal spectral energy distribution (SED) of the HD 181327 disk (See Fig. 4 and Table 2). 
With these data we find a slightly larger thermal IR excess, L^/L^ = 0.25%, than Manning 
& Barlow (1998), and a somewhat warmer blackbody equilibrium temperature for the disk 
grains than the 65 K suggested by Zuckerman & Song (2004). A least-squares blackbody 
fit to the 24 /im and longer flux densities indicates an equilibrium temperature of 77 K ± 5 
K for the thermally emissive grains (after subtracting off the small underlying contribution 
from the stellar photosphere). Following Zuckerman & Song (2004), the warmer grain 
temperature would also suggest - for grains emitting as blackbodies - a smaller characteristic 
orbital radius of « 22 AU. The 12 /zm emission detected by IRAS is not in excess of the 
stellar photospheric level (see Fig 4) as confirmed from a recent Spitzer/IRS spectrum of HD 
181327 (Chen et al. 2006). 

5. THE HD 181327 RING-LIKE DISK 

The NICMOS 1.1 fim PSF-subtracted image of the HD 181327 disk probes closer to 
star than the ACS 0.6 /zm imagery at all azimuth angles (9; with 9 = 0° at the position angle 
closest to the line-of-sight toward the Earth; see §5.1). Because the effective inner radius 
of the ACS image was beyond the peak of the radial surface brightness (SB) distribution 
at most azimuth angles (as revealed by the NICMOS image), we determined the systemic 
geometrical and over-all photometric properties of the HD 181327 disk from the NICMOS 
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image alone. 

The HD 181327 debris disk seen in the NICMOS 1.1 p image is in the form of an 
inclined ring (figure 2A). The radial SB of the ring declines monotonically from the (de- 
projected) radius of its peak brightness, R pea k, both inward (appearing cleared toward the 
center) and outward (appearing bounded or truncated) from R pea k- 

To establish the geometry of the distribution of debris-scattered light in the HD 181327 
ring, we fit image isophotes to nested ellipses and simultaneously solved for the ring ax- 
ial lengths, major axis orientation, and photocenter, by iterative x 2 minimization to the 
NICMOS image. The location of the star (occulted by the coronagraph) was determined 
using acquisition offset slew vectors downlinked in the engineering telemetry applied to star's 
position measured from the well-exposed and critically-sampled, pre-slew target acquisition 
images with a precision of < 1/20 pixel (~ 4 mas; 0.2 AU at 50 pc). 

5.1. Geometry of the Ring 

RING GEOMETRY. Assuming intrinsic circular symmetry, we found the ring inclina- 
tion, i = 31?7 ± 1?6 from a face on viewing geometry, with the disk major axis oriented at a 
PA of 107° ± 2° east of north. The best elliptical fits find the radial peak in surface bright- 
ness at a semimajor axis length of 22.5 ± 0.3 pixels (1"706 ± 0"023). Given HD 181327's 
annual parallax of 19.77 ± 0.81 mas, we find the radius of the debris ring, R pea k = 86.3 ± 
3.7 AU. 

BILATERAL SYMMETRY. Figure 5 shows the disk deprojected from the plane of the 
sky with the major axis on the image horizontal and the location of the occulted star marked. 
The HD 181327 ring surface brightness appears to be bilaterally symmetric, and the disk 
surface brightness peaks, radially, at the same distance at all azimuths. The absence of 
any significant departure from circular symmetry is illustrated in a low-pass filtered image 
(Figure 5B) and the corresponding isophotal map of the deprojected disk (Figure 5C). The 
ring appears to be centered on the star to within a measurement uncertainty of ~ 0.25 pixels 
(1 AU at 50 pc). 

5.2. Photometry of the Disk 

1.1 fim FLUX DENSITY. Because of residual "noise" in the PSF subtraction close to 
the star, there is little useful information at r < 1'.'2 from the star. The flux density from 
the ring dominates over the subtraction residuals at larger radii. We measured the total 1.1 
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/im flux density of the disk at l'/2 < r < 5'.'0 as 9.6 ± 0.8 mJy. The main contributor to 
the ~ 8% measurement error arises from imperfect removal of the underlying stellar PSF by 
under- or over- subtraction due to uncertainties in flux scaling. We assess the magnitude of 
this error from null PSF-PSF subtractions of the reference stars, scaled to the flux density 
of HD 181327 (e.g., see Schneider et al 2001 and SSH05). There may be more dust-scattered 
light interior to this radius. However, as will be shown, the radial surface brightness of the 
disk declines steeply inward from the brightest part of the ring at 1'.'71 (deprojected). The 
small area interior to I" 2 is likely to only weakly contribute to the total light-scattered disk 
flux density. 

1.1 urn SCATTERING FRACTION. Our direct (non-coronagraphic) images of HD 
181327, were saturated in their PSF cores (as expected). Hence, we transform the 2MASS 
J=6.20 magnitude of HD 181327 to the NICMOS F110W band using the CALCPHOT task 
in the STSDAS/SYNPHOT 5 synthetic photometry package. Using an F6V input spectrum 
we find the stellar flux density in NICMOS F110W band is 5.9 Jy, hence, the fraction of 
starlight scattered by the disk, f scat , at r > l'/2 is 0.170% ± 0.015% of the total stellar 
light. We include in our error estimation the uncertainty in our flux density measure, and 
the absolute photometric calibration and transfer of the stellar flux density from 2MASS/J 
to NICMOS/F110W. This 1.1 /im scattering fraction is roughly comparable to the 0.25% 
thermal IR excess found from the > 24 /im color-corrected IRAS and Spitzer flux densities 
(§4). 

1.1 AXIAL RADIAL SURFACE BRIGHTNESS. The SB was measured in radial 
"strips" along the disk major and minor axes, as well as in the direction along the image 
Y axis in the orientation of the ACS/HRC image (see Figure 6) with radial increments of 
75.8 mas (one NICMOS camera 2 pixel) in both outward directions from the occulted star. 
At each radial increment the SB was determined from the median flux density measured 
in regions 1 pixel (along the radial) by 0'/5 (6.5 pixels, corresponding to a ± 10° sector 
perpendicular to the radial at R pea k)- The major axis SB profile is symmetrical about the 
star with a peak SB of ~ 0.97 mJy arcsec" 2 at r = l'/71. In comparison, the asymmetrical 
SB peaks along the minor axis on opposite sides of the star (~ 1.37 and 0.79 mJy arcsec -2 ) 
are strongly indicative of directionally preferential scattering along the line-of-sight. 

1.1 DEPROJECTED A ZIM U TEA L - MEDIA N RADIAL SB. We measured the disk's 
1.1 fim SB profile as a function of distance from the star with elliptical aperture photometry. 
In doing so we assumed intrinsic circular symmetry and a disk inclination of 31?7 from a 
"face-on" geometry. The measurement apertures were elliptical annuli of the same aspect 



5 http://www. stsci.edu/resources/software_hardware/stsdas/synphot 
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ratio of the apparent disk ellipse (1:0.851) and co-aligned with the disk axes. The elliptical 
annuli were 75.8 mas wide on the major axis, and were incremented in semimajor axis length 
by the same amount. The deprojected radial SB profile medianed over all azimuth angles 
is shown in Figure 7. The radial SB was found to peak at all azimuths at isophotal radii 
of 86.3 AU from the star (i.e., 1'.'706 if deprojected). The SB of the ring medianed over all 
azimuths at this peak radius was found to be 1.00 ± 0.07 mJy (F110W ma9 = 15.62 ± 0.07) 
arcsec -2 (see figure 7), very closely corresponding to the ansal brightness of the ring on the 
major axis. 

RING WIDTH. On the azimuthal median, the measured FWHM of the deprojected ring 
(see figure 7) is / /538 (~ 29 AU). At all azimuth angles the SB of the ring declines from the 
1'.'706 (deprojected) radius of peak brightness more steeply toward the interior. The e -1 SB 
falloff from the azimuthally medianed peak is 0"334 (16.9 AU) inwardly and 0"380 (19.2 AU) 
outwardly. In the region l'/2 < r < 5'.'0 approximately 70% of the 9.6 ± 0.8 mJy scattered 
light flux density at 1.1 /im is contained within annulus ± 18 AU in width centered on the 
brightest zone at r = 86 AU from the star. 

0.6 fim SURFACE BRIGHTNESS. Because of the strong azimuthal and radial depen- 
dencies in the ACS PSF-subtraction residuals, SB profiles in the region of the ring were 
measured from the ACS/HRC 0.6 /im image only at r > l'/4 in azimuthal sectors within ± 
10° of the image Y axis. For reasons discussed in §3.2, these are the only sectors in this an- 
nular zone where the flux density of the disk could be measured with reliable error estimates 
against the azimuthally anisotropic "background" of diffractive PSF residuals. To permit a 
direct comparison to the NICMOS 1.1 /jm image, the ACS image was resampled to the NIC- 
MOS camera 2 pixel scale. The SB was measured in both directions from the occulted star 
along the ACS/HRC image Y axis in manner identical to the equivalent NICMOS measures. 
The ACS/HRC image Y axis intersected the HD 181327 ring at 22° from the disk major 
axis as determined from the NICMOS image (see figure 6). The 0.6 /im radial SB profile 
peaks at ~ 1"74 (at P.A. = 95°) and l'/66 (at PA = 275°), with a mean SB of 0.57 ± <0.17 
mJy (F606W ma9 = 16.9 ± < 0.3) arcsec -2 , declining in intensity more rapidly inward than 
outward. 



5.3. Photometric Model of the 1.1/xm Ring-Like Disk 

Without consideration of the physical properties of the disk grains, we constructed a 
photometric model of the disk, predicated on the NICMOS 1.1 fim image, to determine 
the width of the ring, the grain scattering efficiency as a function of azimuth angle (i.e., 
its scattering "phase function" in projection), and to test the measurement of its isophotal 



11 



radius. 



FORM OF THE MODEL. We built a photometric model of the disk assuming an 
inclined, intrinsically circularly symmetric, ring 6 . Based on the results of the isophotal ellipse 
fitting (§5.1), we fixed the disk inclination to 31?7, and orientation to 107°. Allowing for the 
r~ 2 dilution of the central starlight, the surface brightness of the ring was then parameterized 
by a set of four free parameters: guc-, Rpeak, SB(#=0°, R pea k), and W, as follows: 

AZIMUTHAL DEPENDENCE. We assumed a sky-plane projected azimuthal angle (0, 
measured in the plane of the disk) dependence in the scattering efficiency (gHc) of the disk 
grains with the line-of-sight to scattering-phase illumination angle, (p = cos _1 (sini cos 9). 
The spatially resolved deprojected surface brightness, SB(#, r), with r measured in the plane 
of the inclined disk, is represented by a scattering phase function of the type suggested by 
Henyey & Greenstein (1941): 



For the inclination of the HD 181327 disk, scattering phase angles in the range of 58.3° 
to 121.7° are explored. 

RADIAL DEPENDENCE. We modeled a radial SB dependence with inward and out- 
ward exponential decays in SB from R pea k- We characterize the ring width, W, as the distance 
between the radial points where the SB at any azimuth, SB(#, r) declines to SB(9, R pea k) x 
e -1 , i.e., 36.8% of the peak SB at the same azimuth. 

CONVOLUTION. The disk models were convolved with the NICMOS camera 2 F110W 
point-source PSF for the SED of an F6 star using a noiseless, high fidelity, model PSF 
generated with the HST "TinyTim" optical model 7 (version 6.3) as described by Krist & 
Hook (1997). 

FITTING PROCESS AND REGION. The set of parameters which best fit the obser- 



6 Wc tested other functional forms, including model disks with uniform and radial power law surface 
brightness dependencies, but only the ring model provided a meaningful fit to the observational data. Both 
the uniform and radial power law SB models "predicted" excessive levels of scattered light interior to the peak 
radius of the ring that was recovered with equal fidelity after model implantation into the target-minus-PSF 
and PSF-minus-PSF images. I.e., if the HD 181327 scattered light ring exhibited such a form, dust-scattered 
light interior to ring in excess of the declining levels actually detected (eg., see fig. 7) would have been 
observable (and to even smaller radii due to more favorable image contrast) with statistical significance. 

7 http://www. stsci.edu/software/tinytim/tinytim. html 
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vation for % = 31?7, PA = 107°, was found by iterative x 2 minimization in observed minus 
model subtractions in radial zone l'/2 < r < 3'.'0. 

The PSF subtraction residuals in the final image have complex spatial characteris- 
tics which increase in intensity with decreasing angular distance from the coronagraphically 
occulted star. Their amplitudes set the inner distance into which these observations can 
effectively probe. In the fitting process we exclude data at r < 1'.'2, as this region conveys 
little (or no) meaningful information. 

At angular distances > 1"2, the instrumental^ induced PSF-subtraction residuals are 
typically about a pixel in azimuthal extent but are correlated on much larger spatial scales 
radially. These PSF "tendrils" impede the visibility of azimuthal variations in the ring surface 
brightness on comparable spatial scales. To illustrate the underlying form of the ring, we 
removed these PSF subtraction residuals by applying an elliptical low-pass filter to the image 
of the ring (figure 2C). We used a ± 9° rotational "boxcar" kernel of 3.5 pixels width (i.e., 
2.5 resolution elements at the radius suggested by isophote fitting), so only structures with 
approximately half that width are passed. Any "real" features of comparable azimuthal 
spatial frequencies are indistinguishable from PSF subtraction artifacts and also smoothed 
in this process. These high angular spatial frequency features are removed from the filtered 
disk image at the expense of a degradation in azimuthal spatial resolution. The two bright 
spots in the NE sector of the ring (displayed as "hard white" in figure 2C) coincide with 
the location of the wider HST+NICMOS diffraction spikes in two of the four original (pre- 
combined) PSF subtracted images. These regions of enhanced brightness likely arise from 
larger spatial scale, but localized, imperfections in the PSF subtractions and are not intrinsic 
to the disk. For this reason, in the fitting process we also exclude these two regions in the NE 
quadrant of the ring due to the likelihood of their arising from large spatial scale artifacts 
in the PSF-subtracted image. 



5.4. Model Fitting Results 

PEAK RADIUS. R pea k was found to be 22.5 ± 0.4 pixels, in agreement with the radius 
found from fitting ring isophotes. This corresponds to a physical radius of 86.3 ± 3.9 AU. 

CHARACTERISTIC WIDTH. W was found to be 8.26 ± 0.9 pixels. In this charac- 
terization, 69.2% of the light scattered by the dust originates in a 35.8 AU ± 4.0 AU wide 
annulus centered 86.3 AU ± 3.9 AU from the star (The 4.1% uncertainty in the distance 
to HD 181327 from its parallax measure has been added in quadrature to the formal fitting 
error in expressing the uncertainty in R pea k and W in AU) . 
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The W = 35.8 AU ± 4.0 AU fit to the ring width (as characterized) is in very good 
agreement with the estimate derived from the direct measures of the median radial profile of 
the deprojected disk. Neither change significantly with azimuth when measured in 9° wide 
azimuthal sectors (in the absence of subtraction artifacts) around the ring. The large scale 
azimuthal uniformity of the FWHM r ; n9 suggests a constant product of the optical depth and 
grain albedo (r x oo) with azimuth, and perhaps, then, azimuthally undifferentiated grains. 

SCATTERING PHASE FUNCTION. The HD 181327 ring exhibits strong directionally 
preferential scattering with maxima and minima along the line of sight, i.e., in the direction 
of the ring minor axis, with gnc — 0.30 ± 0.03 as found in the fitting process. (Discrete SB 
measurements along small azimuthal sectors at constant deprojected radii depart from this 
fit due to biasing by local PSF-subtraction residuals 8 ). A photometric model of the disk, 
convolved with the NICMOS camera 2 PSF, g HG = 0.30, and i = 31?7 from face-on (shown 
in Figure 2B) predicts ansal SBs of 0.863 mJy arcsec -2 and minor axis SB peaks of 1.37 and 
0.56 mJy arcsec -2 , each in agreement with the photometric measures (figure 6) within the 
1 a measurement uncertainties assessed w.r.t. the azimuthal median. 



6. NICMOS + ACS 

The NICMOS and ACS instruments each better probed different spatial regimes of 
the HD 181327 circumstellar disk (for a comparative overview of the HST coronagraphs 
see Schneider et al. (2003) 9 ). While the spatial domain of overlap with high photometric 
fidelity in the radial region of the debris ring is small, for the reasons noted in §3.2, together 
the NICMOS and ACS observations further inform on the global structure and properties 
of the disk, and its environment. The NICMOS imagery reveals the detailed structure of 
the disk and region interior to ~ 3", whereas the larger ACS field and deeper integrations 
reveal the faint outer nebulosity surrounding the ring at r > 4". In the regions commonly 
sampled the nearly factor of 2 wavelength coverage afforded by the ACS/F606W (O.Qfim) and 



8 Thc SB of the disk at R pea k was also discretely measured as a function of elliptical azimuth around the 
ring. Measurements were made in "face on projection", with the disk surface brightness preserved in the 
geometrical transformation. Circular photometric apertures of both 3.5 and 7 pixel diameters (one and two 
times the FWHM of the SB peak) were used to measure the SB every 9° (40 measures), and every 18° (20 
measures) around the ring, respectively. Aperture dilution of the non- uniform flux density distribution is less 
severe in discrete measurements with the smaller aperture, but these individual measures are more influenced 
by local subtraction residuals than with the larger aperture. The SB measures using both apertures were in 
agreement with a gHG = 0.30 ±0.03 scattering phase function. 

9 http://nicmosis.as.arizona.cdu:8000/POSTERS/AAS JAN2004.CORON.jpg 
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NICMOS/F110W (1.1/xm) filters provide rudimentary color information on the disk grains. 

6.1. [F606W] - [F110W] Color of the Disk at R peak 

HD 181327 has a V-J (similar to [F606W]-[F110W] for an F5/F6V star) color index of 
~ +0.80. As noted in §5.2, the 0.6 fxm SB of the disk could be measured reliably only in 
small regions along the debris ring between 174 and 2'/0 flanking the ACS/HRC image Y 
axis (indicated in figure 3A) at celestial position angles of 95° and 275°, 22° from the disk 
major axis (as shown in figure 6). Identical measures from the best-fit 1.1 /zm scattered light 
model of the ring based upon the NICMOS image found SBs of 1.30 and 0.72 mJy arcsec -2 
at the same position angles, respectively. The mean 1.1 /im SB in these regions corresponds 
to 15.61 ± 0.08 (F110W ma9 arcsec -2 ). Hence, in these regions where the disk-scattered light 
peaks radially, the disk grains appear significantly red with a [F606W] - [F110W] color index 
of +1.32 ± 0.3. Correcting for the color of the star, the disk grains in this region have an 
intrinsic [F606W] - [F110W] color index of +0.5 ± 0.3. 

6.2. The Outer Regions of the Disk 

Beyond R pea k, in the region of high photometric reliability of 100 - 200 AU (r < 4" in the 
deprojected image NICMOS image), the 1.1 /im SB declines as 15.75 x r -5 - 2 mJy arcsec -2 
on the azimuthal median (see fig 8). Using the ring inclination and orientation derived from 
the NICMOS image, we similarly deproject the ACS image and find a very similar power 
law index with a monotonic decline in the 0.6 /im SB of 7.98 x r -4 9 mJy arcsec -2 in the 
same region (also shown in figure 8). The azimuthally medianed [F606W] - [F110W] color 
index of the disk throughout the 100 - 200 AU region is statistically identical to that found 
at the two measurable locations at 86 AU (§5.2). Hence, from these limited data we see no 
evidence of radial differentiation in the color of the disk grains. 

The ACS image reveals a diffuse low-SB nebulosity to a distance of < 9" from the 
star (figure 9, left panel), the inner portion (beyond ~ 4") of which is only hinted at in 
the NICMOS image (figure 9, right panel). The nebulosity seen in the ACS image, with a 
mean 0.6 /im SB in the 4" < r < 9" region of ~ 21.5 (F606W ma9 ) arcsec -2 , is azimuthally 
asymmetric with the brightest part in an ~ 90° sector to the North of the star and "curving" 
in the outer regions toward the east. HST coronagraphic observations have revealed spiral 
structures in the circumstellar disks of the Herbig AeBe stars HD 100546 (Grady et al. 2001) 
and HD 141569A (Clampin et al. 2003). In the latter case, Clampin et al. postulated that 
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HD 141569's outer "arms" might have arisen from tidal interactions from a passing star, 
which could also be the causal mechanism for both the disk's stellar photocentric offset 
and its large azimuthal surface density variations. Neither of these phenomena are seen in 
the HD 181327 disk. The southwest arm of the HD 141569A disk seems to be associated 
with its M-star companions, but no tidally interacting companions have been found for 
HD 181327 (see § 6.3). Kalas, et al. (2002) have reported on the existence of Pleiades-like 
reflection nebulosities in the environments of Vega-like stars of spectral types K2V-B9.5V 
(i.e., similar to HD 181327), though an examination of Spitzer 24 and 70 /mi images of the 
field around HD 181327 (§ 4) does not reveal any emission structure the could be identified 
with the low-SB feature seen in the ACS 0.6 fxm image. Hence, the observations we report 
here, unfortunately, do not inform on the physical nature of the axially asymmetric nebular 
structure in proximity to HD 181327 that may, or may not, be physically associated with its 
disk. 



6.3. Tests for Companionship 

In §1 we suggest that HD 181327 and the 10-30 Myr star HR 7329 (Lowrance et al. 
2000) are a common proper motion pair, strengthening the argument for HD181327's youth. 
Four point-like objects (stars) are seen in the more immediate vicinity of HD 181327 in 
both the NICMOS and ACS images (denoted a-d in Figure 9). Although stars a and b 
were within the 1276 x 1276 SHARP-I (Neuhauser et al. 2003) and 1277 x 1277 SHARP-II 
(Chauvin et al. 2003) survey fields, neither were detected in either study. We also searched for 
low-luminosity (sub-stellar) companion candidates in close angular proximity to HD 181327 
using the NICMOS coronagraphic images. Following Schneider & Silverstone (2003) and 
Lowrance et al. (2005), we implanted model "PSF stars" (using noiseless TinyTim 6.3 PSFs 
(Krist & Hook 1997)) into a difference image made from of the two differentially-oriented 
coronagraphic images of HD 181327. Even in the presence of the HD 181327 disk, at and 
interior to the radius of peak surface brightness, we found (similar to Weinberger et al. 
(1999) in the case of HD 141569 A) that at this approximate age, brown dwarf or giant 
planet companions with masses > 10 Mj upiter would be detectable in to an interior radius 
of pa 1". To these approximate sensitivity limits, no close-proximity sub-stellar companion 
candidates were found. 

The age of HD 181327, presumed ~ 12 Myr as a member of the f3 Pictoris moving group, 
could be further constrained by estimating the age of a low- mass (stellar) companion (e.g., 
by the presence of Li, H Q equivalent width, X-ray activity, etc.). We tested the possibility 
of companionship of the four stellar-like objects appearing in the NICMOS and ACS images 



-16- 



through differential proper motion (PM) measures over the 1.5 years between observations. 
On the aggregate, the four had measured differential PMs of (ARA, ADEC) = (-42.3, +83.9) 
mas yr -1 with a 1 tr dispersion of ± (3.3, 16.0) mas yr~ x w.r.t. HD 181327, suggestively 
indicative of the reflective motion of HD 181327 (the largest "outlier", star c, was discrepant 
from the ensemble mean by (4.8, 23.6) mas yr _1 ). The proper motion of HD 181327 itself 
is given as (+24, -82) mas yr _1 by Perryman et al. (1997). Given the very good internal 
agreement in the measured PM of the four stars, the dispersion being < (0.14, 0.20) a of the 
PM of HD 181327, all four - unfortunately - are very likely background stars 10 . 



7. DISK/GRAIN PROPERTIES 

7.1. Inversion of the Scattered-Light SB Profile: Approach 

We modeled the scattered light surface brightness, SB(x), of an axisymmetric, inclined, 
and optically thin disk, measured at a projected distance x from the star along the apparent 
major axis, following the methodology as described by Augereau & Beust (2006) such that: 



<P* J x r l — x 2 cos(z) y/r 2 — x 2 

where i is the disk inclination from pole-on, p{r, z) is the disk number density in the 
cylindrical frame centered onto the star and attached to the disk, r max is the outer disk 
radius and $* is stellar flux at Earth at the considered wavelength. a sca and /(</?) are the 
scattering cross section of the grains and the anisotropic scattering phase function at phase 
angle ip, respectively, both averaged over the grain size distribution. 

In the case of HD 181327, the mean phase function /(</?) is observationally constrained 
and can be approximated by a Henyey-Greenstein type phase function with guc = 0.3 ± 



10 The magnitude of the apparent discrepancy in the reflective PM in RA of HD 181327 over the 1.5 
years between the observations, assessed from the four star ensemble compared to its a priori measured 
value, is ~ 28 mas (approximately 1/3 of a NICMOS camera 2 pixel). We attribute this discrepancy to an 
(unresolved) inter-camera differential astrometric systematic uncertainty, likely due to differential camera 
distortion corrections and/or occulted target centration errors (as previously noted in §3.2 for the ACS 
observation). We note that all four "reference stars" were spatially distributed within a single quadrant (to 
the west of HD 181327) due, in part, to the intersecting areas of the two camera fields. A more isotropic 
distribution of reference stars (or a two-epoch observation with the same instrument) would have provided 
a better absolute astrometric solution. A fifth, brighter, star appears in the ACS field (to the east of 
HD 181327) only, so its differential PM could not be measured. 
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0.03 at 1.1 jum, as noted in §5.4. We also presume an axisymmetric disk predicated on the 
observed bilateral symmetry of the scattered-light ring in the NICMOS image. Therefore, for 
a given vertical dust distribution, Z(r,z), the product of the midplane number density, p(r, 
0) = p[r,z)j Z(r,z), and the mean scattering cross section, a sca , can be obtained by directly 
inverting the observed surface brightness profile SB(x). We do so using the technique detailed 
in Augereau & Beust (2006) and discuss the properties of the dust ring as a function of 



We characterize the vertical distribution with a three parameter function: Z(r,z) = exp(- 
(\z\/ H(r)y'); where H = ifo(?Y r o) is the vertical scale-height and r$ is the observed radius 
of the ring (86.3 AU). We calculated the midplane number density for both exponential (7 
= 1) and Gaussian (7 = 2) profiles assuming either linearly flaring (f3 = 1) or a nearly flat 
((3 = 0.1) disks. In each case, twenty scale-heights, 1 AU < H < 20 AU were considered 
corresponding to 0.012 < H /r < 0.23 (see §7.3). 

To compare the profiles, we introduce the disk surface density, £(r), that can be written 



where r±(r) is the vertical optical thickness at 1.1 /im and u is the mean albedo in the 
NICMOS F110W passband averaged over the grain size distribution. C 7 is a constant value 
depending upon 7 (C\ = 2, and C2 = s/k). 



Figure 10 shows several a sca £(r) = u r±(r) profiles obtained by directly inverting the 
averaged NICMOS F110W SB profile. For sufficiently large values of Hq, the vertical optical 
thickness at 1.1 /im is independent of the shape of Z(r,z) while the width of the ring 11 of 
43 ± 1 AU (FWHM ~ 32.5 ± 0.5 AU) depends slightly on the vertical profile. At large 
distances, the surface density profile falls off as ~ r -3 , in agreement with the observed r~ 5 
SB profile. A surface density profile E(r) oc r -3 is fully consistent with dust grains supplied 
by parent bodies on circular orbits and placed in orbits of high eccentricity by a drag force 
oc a 2 / (r 2 + z 2 ) where a is the grain size (e.g., see Eq. 3 and Fig 1 of Lecavelier des Etangs et 
al. (1996)). This model additionally requires a differential grain size distribution that does 
not depart too much from a a~ 3 ' 5 power law, and is in good agreement with the conclusions 



full width at c 1 of the maximum 



Z(r,z). 




a sca p(r, z) dz = C 7 H(r) a sca p(r, 0) = u T±(r) (3) 



7.2. Inversion of the Scattered-Light SB Profile: Results 
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from the SED fitting (see §8). Moreover, because the model needs the dust dynamics to be 
largely unaffected by gas, the shape of the observed dust profile may be indicative of low gas 
density as it has been similarly shown for (3 Pictoris (Thebault & Augereau 2005). 

As HD 181327 is a star of spectral type F5/F6V, radiation pressure contributes, at least 
partly, to the total drag force on the grains. In that picture, the dust particles observed 
at distances larger than ~ 100 AU most likely originate from a belt of large solid bodies 
peaking close to 87.6 AU. 

The midplane optical thickness, ry, in the NICMOS F110W passband is related to the 
vertical thickness through 



w T|1 (r) = [ asca p{R ' 0) dR = [ 70§) dR (4) 

with T|| being the 1.1 jum optical thickness calculated from the region between the 
inner edge of the disk to a distance r. As shown in Figure 11, with H too small the disk 
appears marginally optically thick in the midplane. Using the constraints on the grains for 
astronomical silicates (Weingartner & Draine (2001); §8), we estimate the mean albedo, u 
to be ~ 0.6. Although the albedo is not well constrained, we can safely conclude that Hq 
must be > 4 AU {Hq/tq > 0.05) to ensure the disk is optically thin in all directions at 1.1 
fim. When this condition is fulfilled, the ring vertical thickness, T±(r), reaches a maximum 
of ~ 0.01 (assuming u ps 0.6) at 87.6 AU (Fig 10). Assuming a grain mass density of 3.5 
g cm -3 and a a -3 - 7 differential grain size distribution between 1 /iin and 1 mm, we infer a 
total dust mass of ss 4.1 M Moon (ps 0.05 M e ). 

The collision time-scale may be estimated as suggested by Backman & Paresce (1993), 

, °~sca f / - \ 

trail — n X 1 5 

vr < a 2 > 2 uo r ± (r) VGM 1 ' 

where ir <a 2 > is the mean geometric cross-section of the grains averaged over the grain 
size distribution, G is the gravitational constant, and M is the stellar mass (~ 1.35 Mq for 
an F5/F6V star). The size distribution derived in § 8 gives o- sca /n <a 2 > ~ 1.4. For such 
grains located within 150 AU of the star, the collision timescale is < 10 5 yr, except very close 
to the inner ring edge where the density is extremely low. At the position of the peak surface 
density, t co u 10 4 yr (hence about three orders of magnitude smaller than the star's age). 
The ring is then collision-dominated and the smallest of the released grains can populate the 
outer disk thanks to radiation pressure commensurate with the r -3 surface density profile 
at large distances. The short collision time-scale implies that any dynamical process for 
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for the dust grains lasting longer than 10 4 " 5 yr can thus be neglected. This is the case of 
Poynting- Robertson drag for which the depletion timescale at 87.6 AU for just-bound grains 
compares to the stellar age (i.e., ~ 10 Myr). I.e., dust cannot populate the disk interior to 
the ring due to Poynting-Robertson drag. 

7.3. Scattered-Light Models and Fits 

Synthetic scattered-light images of the HD 181327 circumstellar disk were generated 
employing the three dimensional radiative transfer code, MCFOST, described by Pinte et 
al. (2006). The product of density, p(r,z), and scattering cross section, cr sca , resulting from 
the inversion of the surface brightness profiles (§§7.1 & 7.2) were used assuming a Henyey- 
Greenstein scattering phase function with gnc — 0.3 (§5.4). Photon packages escaping 
the system with inclination angles between 30?5 and 32?5 were recorded in the synthetic 
images. Representative model disks for a subset of characterizing parameters a,/3, and H , 
corresponding to several of the profiles discussed in §7.2 (and presented in figures 10 and 
11), are shown in figure 12. 

The spatially resolved flux densities predicted by the MCFOST scattered-light models 
were fit to the NICMOS image in the annulus of high SNR (r l'/3 — 2'/3) enclosing the 
observed scattered-light ring. Each of the previously discussed models, tiling the parameter 
spaces illustrated in Figure 12, were fit to the NICMOS image by \ 2 minimization, scaling 
the intensities of the predicted scattered-light patterns to the observed spatially-resolved flux 
densities. A representative example (the /3 — 1, 7 = 2, H — 1AU case), and its fit to the 
NICMOS image after \ 2 minimization, is shown in Figure 13. Within each 7) "family" 
of models, the reduced x 2 estimation was relatively flat for H < 8 AU. In §7.2 we separately 
concluded that H must be > 4 AU (H /r > 0.05), hence we suggest 4 AU < H < 8 AU. 

8. GRAIN SIZE DISTRIBUTION 

Both the ACS (0.6 /jm) - NICMOS (1.1 /im) color of the HD 181327 disk and the 
measured scattering asymmetry factor (gHG — 0.3 at 1.1 /im) indicate that the scattered 
light images are dominated by sub-micron grains. The grain properties that may explain the 
scattered light observations should also be consistent with other independent constraints for 
the grains. In particular, grains at the observed radius of 86.3 AU should produce an infrared 
signature consistent with the measured long-wavelength excesses and the grains properties 
should be consistent with basic dynamical considerations. Because HD 181327 is an F5-F6 
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type star, the radiation pressure force on the grains is expected to be significant and can 
expel sufficiently small grains from the disk on very short time-scales. 

We have assumed a unity normalized differential grain size distribution ( f amax dn(a) = 

1) of the form dn(a) oc a K da between a minimum grain size, a min , and a maximum grain size, 
Umax- We thus search for (a m i n , «) pairs consistent with the four constraints of color, gnc 
IR excess, and blow-out size limit (thus, corresponding to the depletion timescale), identified 
above. In the following, the maximum grain size will be assumed to be large enough to not 
affect the results. The grain optical properties are calculated with the Mie theory valid for 
hard spheres and for astronomical silicates (optical constants from Weingartner & Draine 
(2001)). 

The measured mean scattering asymmetry factor, g HG = 0.3 ± 0.03 in the F110W 
passband, then reads 

/Q>max 
g(a) tx a 2 Q sca dn(a) / a sca (6) 
-min 

where Q sca is the dimensionless scattering efficiency, a sca = f amax na 2 Q sca dn(a) and g(o) 
is the grain-size dependent scattering asymmetry factor at 1.1 /im. As the disk is optically 
thin, the measured disk color is dependent upon the ratio of scattering cross sections and 
phase functions at the two observed wavelengths. The mean scattering asymmetry factor 
at 0.6 /im, unfortunately, could not be measured from the ACS image. We, therefore, have 
made the rough assumption that the phase function in the ACS passband is comparable 
to that found from the NICMOS (1.1 /im) observations. In that case, the measured color 
directly corresponds to the mean scattering cross section ratio: [F606W] — [F110W] = 2.5 log 
10 (<7 aco [F606W]/<7 aco [F110W]) = +0.5 ± 0.3. 

To further constrain the grain size we used the surface density profile obtained in §7 to 
calculate infrared flux densities in the A > 12 /im IRAS 12 and Spitzer/MIPS bands (Table 

2) as a functions of a min and k using the Augereau et al. (1999) model. Fitting the SED 
leads to reduced x 2 contours in the (a min , k) parameter space as illustrated in Figure 14. 
The best fit is obtained for minimum grain sizes close to 1 /im, in very good agreement 
with the radiation-pressure induced minimum grain size. The SED fitting also reduces the 
range of possible k values, the best fit being obtained for a power law index of as -3.7, quite 
close to the theoretical -3.5 power law index that holds for systems in collisional equilibrium 
(Dohnanyi 1969) (dashed line in Figure 14). 



12 Thc 12 /im emission detected by IRAS arises from the stellar photosphere, not the disk. 



- 21 - 



The ACS-NICMOS disk color is reproduced with only slightly smaller (0.4 < a min < 
0.6 jum) grains than those with sizes that best fit the thermal SED (1 < a min < 2 jum). 
Nevertheless, these grains remain too large to reproduce the measured asymmetry factor 
Zhg (requiring a m i n < 0.1 /im grains in our formulation). In the frame of the Mie theory, 
the observed mean asymmetry factor can indeed only be reproduced when scattering is 
dominated by very small grains which is achieved when both a m i n and k are extremely 
small. This inconsistency may relate to the fact that while we have calculated the asymmetry 
factor for hard silicate spheres (a simplifying assumption inappropriate for grains with more 
complex structures), the scattering phase function is very sensitive to the shape and structure 
of the grains. Additionally, the HD 181327 disk may possess disk grains that: (a) originate 
in more than just a single population, (b) are radially segregated in distance from the star, 
and (c) have a grain-size distribution not well-represented by a single power-law. 

The grains which are predominantly responsible for HD 181327's thermal infrared excess 
may not be those which give rise to the circumstellar light scattered at 0.6 to 1.1 fira. In 
the case of the debris system circumscribing HR 4796 A (Schneider et al. 1999), Augereau 
et al. (1999) suggested a two-component dust population with large grains of ISM-like 
composition (moderately porosity amorphous silicates) at the radius of the scattered-light 
disk, and high porosity crystalline grains close to the star. The thermal emission from 
HD 141569's circumstellar disk was shown by Marsh et al. (2002) to arise from within 1" 
(100 AU) of the star from 12 - 21 fj,m imaging, while light scattered by the disk at 1.1 and 1.6 
/im (Weinberger et al. (1999); Augereau et al. (1999b)) originates (and extends) significantly 
further out. While HD 181327's thermal SED is not inconsistent with the possibility of ISM 
like grains (e.g., k -3.5 with a min < 1 /im ; Mathis, Rumpl & Nordsieck (1977)) in the 
scattered-light debris ring, the presumption of a single population of hard, spherical, silicate 
grains would be inapplicable for grains which have grown (evolved) from ISM-like structures, 
and may be spatially (radially) segregated in distance from the star. 



9. DISCUSSION 

The inferred spatial distribution of dust detected about stars, from thermal infrared 
excesses (as measured from IRAS and Spitzer), often predicts central clearings or "holes" to 
explain the observed spectral energy distributions. Scattered-light imagery of circumstellar 
debris directly provides constraints on the spatial distribution of dust independent of the dust 
grain temperature in regions sufficiently far from the central star which can be instrumentally 
probed with high spatial resolution. The central regions of dusty circumstellar disks that are 
nearly edge on to the line-of-sight (e.g., (3 Pictoris, AU Microscopii and HD 32297) remain 
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hidden at optical and near-IR wavelengths due to the opacity of the dust. The inner regions 
of debris disk systems which are more inclined to the line-of-sight, such as HD 181327, 
are revealed directly and give support to the interpretation of inner holes from analyses of 
debris disk SEDs. While only a very small number of light-scattered images of circumstellar 
debris disks exist, the growing number allows inter-comparison of the disk characteristics 
with stellar type and age (when sufficiently constrained; see Table 3 and Figure 15). With 
the notable exception of a PsA, which thus far has only been observed in scattered light with 
HST/ ACS, all of the debris disks listed in Table 3 have been observed with scattered-light 
imagery in the optical and the near-infrared. Very recently, two newly discovered, scattered- 
light disks around older (0.3 - 1 Gyr) stars, HD 53143 and HD 139664 (Kalas, et al. 2006), 
have been reported, adding importantly to this still very limited sample. 

Of the debris disks imaged to date, the HR 4796A ring-like disk is, in many respects, 
morphologically most similar to HD 181327's. Both scattered-light disks appear as relatively 
narrow annuli, have similar near-IR scattering fractions and thermal IR excesses, and phys- 
ically are of similar size (70 and 86 AU in radii, respectively). Both stars, while not coeval 
(HR 4796A being a member of the TW Hya association, and HD 181327 in the f3 Pictoris 
moving group) are of similar ages (~ 10 Myr), and spectral types. 

Directionally preferential scattering, symmetrical about the disk minor axis (i.e., with 
respect to the line-of-sight), is seen in the HD 181327 disk. Such non-isotropic scattering 
might be expected from circumstellar particles evolved from a primordial ISM-like population 
(comparable in size to, or larger than, the 1.1 /zm wavelength of the observation), and similar 
behavior is seen in the HR 4796A debris ring, and in a PsA's and HD 107146's as well. 

These four ring-like debris disks vary by factors of ~ 2 in radius (from the already 
noted smaller HR 4796A and HD 181327 to the larger HD 107146 (133 AU; Ardila et al. 
(2004)) and a PsA (145 AU; Kalas, et al. (2005)), and ~ 4 in ring width-to-radius (w:r) 
ratio (from ~ 0.17, for both HR 4796A (Schneider et al. 2004) and a PsA, to 0.64 for 
the very broad HD 107146 ring, and HD 181327 in between at 0.38). The disks about 
HR 4796A and a PsA, both early A stars, are most disparate in radius but most similar 
in w:r ratios, suggesting a lack of simple scalability in ring geometries conflated (if not 
determined) instead by dynamical and kinematic processes in play. For example, in the case 
of HR 4796A, the outer radius of its disk may be truncated by the dynamical influence of its 
M-star companion (e.g., Artymowicz & Lubow (1994)). The sharp inner and more shallow 
outer edges of the illumination-corrected scattered-light radial profiles of HD 181327 (Fig. 
7) and a PsA (figure 3 of Kalas, et al. (2005)) are remarkably similar, suggestive of similar 
clearing and/or confinement mechanisms for the disk grains. 

The photocenter of the a PsA's dust ring is asymmetrically offset from the position of 
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the central star. As suggested by Kalas, et al. (2005), such an asymmetry may implicate 
the existence of one (or more) planetary perturbers. The HR 4796A debris ring exhibits 
an ~ 20% SB asymmetry at the ring ansae at near-IR and optical wavelengths (Schneider 
2004) (also seen in it's mid-IR thermal emission as reported by Telesco et al. (2000)), with 
no photocentric offset in the position of the star relative to the ring ansae detectable in the 
scattered-light observations. Such a non-axisymmetric SB asymmetry about the ring minor 
axis could also arise from a local enhancement (or depletion) in the dust density distribution 
due to perturbations from unseen, co-orbital, planetary- mass bodies (e.g., Ozernoy et al. 
(2000)). Though, alternative mechanisms, dependent upon the amount of residual gas in 
the system, could account for the ring structure itself (e.g, Klahr & Lin (2005); Takeuchi & 
Artymowicz (2001)) without invoking planetary-mass companions. 

The HD 181327 ring morphology is bi-axially symmetric, centered on the star, and 
devoid of any readily apparent azimuthal SB anisotropies aside from those dependent upon 
scattering phase angle. Such a ring structure could arise from gas drag in the absence of 
planets, though the gas-to-dust ratio in the HD 181327 disk is, currently, unknown. 

10. SUMMARY 

We have imaged a circumstellar debris disk about the ~ 12 Myr old F5/F6V star and (3 
Pictoris moving group member HD 181327 in 1.1 /im and 0.6 jum light scattered by the disk 
grains. Our 1.1 /im PSF-subtracted NICMOS coronagraphic observations reveal a bilaterally 
symmetric ring-like structure, 31?7 ± 1?6 from face-on with a major axis PA of 107° ± 2° 
east of north. We estimate the 1.1 /im disk flux density in the region l'/2 < r < 5'.'0 to be 9.6 
± 0.8 mJy (0.17 % of the total stellar light) with ~ 70% of the disk flux density contained 
within a 36 AU wide annulus centered on the brightest zone at r = 86 AU from the star, well 
beyond the blackbody characterized 22 AU radius of the thermally emissive dust detected in 
the IRAS and Spitzer/MYPS bands. The scattered-light disk appears to be centrally cleared 
(to within 1'.'2 of the star, which is the limiting region which can be effectively probed with 
our NICMOS observations) and centered on the star. 

The surface brightness of the ring of debris-scattered light is well represented by a 
Henyey & Greenstein (1941) type scattering phase function with g HG = 0.3. Modeling the 
disk, we find a vertical scale height H = 5-10% of the 87.6 AU radius of the maximum surface 
density at that radius. The observed dust ring explains the thermal IR excess emission above 
the stellar photospheric level measured by IRAS and Spitzer provided that a min > 1 /im with 
a grain size distribution oc a -3 ' 7 , presuming hard, spherical, silicate grains. With these grain 
parameters, and surface density profiles derived by direct inversion of the SB profiles obtained 
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from the NICMOS imagery, we infer a total dust mass of ps 4 M M oon (~ 0.05 M ffi ) in the 
collision-dominated ring. Such grains, however, fail to reproduce the gnc — 0.3 scattered- 
light asymmetry factor, suggestive of a more complex (possibly mixed) population of disk 
grains. To better constrain the possible grain sizes and their radial distributions, in HST 
cycle 15 we will obtain 2 /iin coronagraphic imaging polarimetry observations with NICMOS 
(a new, recently enabled, observing mode; Hines & Schneider (2004)) of the HD 181327 disk 
(GO program 10847). 

Additional circumstellar light is seen from 2"-4" (100-200AU) in both 1.1 jum NICMOS 
and 0.6 /xm ACS imagery declining in SB as ~ r -5 in both bands. The SB profile beyond 
100 AU is consistent with the blow-out time scale for grains in the inner disk by radiation 
pressure, collisional grains with a size distribution as discussed, and low gas pressure in 
the disk. A low SB diffuse outer halo (V 21.5 mag arcsec~ 2 ) is seen in the ACS image 
extending as far as ~ 9" (450 AU) to the north of the star. 

We are grateful to F. J. Low and P. S. Smith of the Spitzer GTO (PID 72) team for 
providing MIPS imaging photometry of HD 181327. We thank Elizabeth Stobie (NICMOS 
Project, University of Arizona) for her significant contributions to this program through her 
development, implementation, enhancement and support of analysis software used through- 
out the course of this study. We appreciate the commentary received from Eric Mamajek, 
particularly in regard to the age estimation for HD 181327 beyond its "guilt by association" as 
a presumed member of the /3 Pictoris moving group, as reflected in our discussion in §§1 and 
6.3 of this paper. This investigation was based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope Science Institute (STScI), which 
is operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated with programs GO/10177 and 
GTO/9987. Support for program GO/10177 was provided by NASA through a grant from 
STScI. S. Wolf was supported by the German Research Foundation (DFG) through the 
Emmy Noether grant WO 857/2-1. J.C.Augereau, C. Pinte and F. Menard acknowledge fi- 
nancial support from the Programme National de Physique Stellaire (PNPS) of CNRS/INSU, 
France. We acknowledge the use of the SIMBAD database. 
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Table 1. NICMOS OBSERVATIONS OF HD 181739 and PSF Refrence Stars 



Visit a Star b 


Orient 


Typc d 


Filtcr d 


Readout Mode d 


ExpTime e 


Exp f 


Datasets g 


67 HD 181327 


37.570 


ACQ 


F180M 


ACQ (ACCUM) 


0.243s 


2 


N8ZU67010 






CORON 


F110W 


STEP32/NSAMP14 


224s 


3 








/" ' 1 \ 1 "\ / \ TV T 

CORON 


F110W 


STEP32/NSAMP10 


96s 


1 








DIRECT 


F110W 


SCAMRR/NSAMP11 


2.03s 


2 




68 HD 181327 


66.970 


ACQ 


F180M 


ACQ (ACCUM) 


0.243s 


2 


N8ZU68010 






CORON 


F110W 


STEP32/NSAMP14 


224s 


3 








CORON 


F110W 


STEP32 /NSAMP10 


96s 


1 








DIRECT 


F110W 


SCAMRR/NSAMP11 


2.23s 


2 




83 HD 4748 


278.970 


ACQ 


F171M 


ACQ (ACCUM) 


0.309s 


2 


N8ZU83010 






CORON 


F110W 


STEP32 /NSAMP12 


160s 


1 








CORON 


F110W 


STEP32/NSAMP14 


224s 


10 








DIRECT 


F110W 


SCAMRR/NSAMP10 


2.03s 


8 




84 HIP 7345 


198.398 


ACQ 


F187N 


ACQ (ACCUM) 


0.763s 


2 


N8ZU84010 






CORON 


F110W 


STEP32 /NSAMP10 


96s 


1 








CORON 


F110W 


STEP32/NSAMP14 


224s 


10 








DIRECT 


F110W 


SCAMRR/NSAMP10 


2.03s 


8 





a ffi>T/GO 10177 Visit ID number 

b 2MASS mags from Cutri et al. (2003), Spectral types from Perryman et al. (1997): 
[HD 181327; VT = 7.095 ± 0.007, J = 6.20 ± 0.02, H = 5.98 ± 0.04, K = 5.91 ± 0.04, F5/F6V], 
[HR 4748; J = 5.52 ± 0.02, H = 5.58 ± 0.03, K = 5.54 ± 0.02, B8V], 
[HIP 7345; J = 5.49 ± 0.02, H = 5.53 ± 0.02, K = 5.46 ± 0.02, A1V]. 

c Position angle (degrees) of image +Y axis (East of North). 

d See Noll et al. (2004) for NICMOS image/readout modes and filters. 

e Exposure time for each exposure (before combining). 

' Number of exposures 

S HST archive data set ID 



Table 2. Spitzer/MIPS & IRAS Thermal IR Flux Denisties 



Wavelength (^m) a 


Source 


Flux Density (Jy) b 


Uncertainty 


Color-Corrected (Jy) 


12 


IRAS 


0.164 


10% 


0.114 ± 0.011 


23.68 


MIPS 


0.220 


4% 


0.223 ± 0.009 


25 


IRAS 


0.248 


8% 


0.306 ± 0.024 


60 


IRAS 


1.86 


4% 


1.96 ± 0.08 


71.42 


MIPS 


1.58 


7% 


1.73 ± 0.12 


100 


IRAS 


1.72 


12% 


1.70 ± 0.20 


155.9 


MIPS 


0.750 


12% 


0.766 ± 0.090 



a Weighted mean wavelengths for MIPS; see http://ssc.spitzer.caltcch.edu/mips/dh/ 
b IRAS flux densities from the FSC 

c for 70K blackbody grains in all bands except IRAS 12 /im: for 6450K photosphere. Sec: 
http://irsa.ipac.caltcch.edu/IRASdocs/exp.sup/ch6/tabsupC6.htm 
http://ssc.spitzer.caltech.edU/mips/dh/mipsdatahandbook3.2.pdf 
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Table 3. Circumstellar Debris Systems Imaged in Scattered Light 



Star Spec Age (Myr) r(AU) a / n i r (r>0'.'3) L ir /L* Assymetry Image 



AU Mic b 


MO 


12-20 


210 




0.0004 


Wing Tilt 


Fig 15a 


P Pic c 


A5 


12-20 


1600 


0.003 


0.0015 


many 


Fig 15b 


HD 32297 d 


AO 


<10 ? 


>400 


0.0033 


0.0027 


SB/W. Tilt 


Fig 15c 


HR 4796A e 


AO 


8 


70 


0.0024 


0.005 


Azim. 


Fig 15d 


HD 181327 


F5/6 


12 


86 


0.0017 


0.0025 


No, g«0.3 


Fig 15e 


HD 141569A^ 


HAeBe/B9.5 


5 


400 


0.0025 


0.0084 


Yes 


Fig 15f 


HD 1071469 


G2 


30-250 


170 


1.4xlO~ 4 


0.001 


No, gw0.2 


Fig 15g 


a PsA h 


A3 


200 


140 


10~ 6 


5xl0~ 5 


Offset, g;s0.3 


Fig 15h 



a rAu reported is scattered-light maximum extent detected for near edge-on disks (and HD 141569A), radius 
of peak surface brightness for ring-like disks (HR 4796A, HD 181327, and a PsA). 

b AU Mic: Kalas, ct al. (2004). 

c (3 Pic: Kalas & Jewitt (1995), / nir estimated from Kalas, ct al. (2000), L ir /L, from Dccin et al. (2003). 
d HD 32297: Schneider ct al. (2005). 
e HR 4796A: Schneider et al. (1999). 

f HD 141569A: Weinberger ct al. (1999), anisotropics see Mouillct ct al. (2001) & Clampin ct al. (2003). 

g HD 107146: Ardila ct al. (2004), g at 0.8 (im. 

h a PsA: Kalas, et al. (2005), L ir /L* from Dccin ct al. (2003). 
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Fig. 1. — Fig. 1 - The robustness of NICMOS PSF-subtracted coronagraphy in revealing the 
HD 181327 circumstellar debris ring is illustrated by these individual (A-D) representative 
astrometrically registered, and flux-scaled difference images, using two different PSF tem- 
plate stars (HIP 7345 in panels A & B; HIP 4748 in panels C & D) with HD 181327 observed 
at two field orientations differing by 29?9 (absolute orientations w.r.t. North as indicated). 
Panel E shows a same-scaled and registered PSF-minus-PSF subtraction of the two template 
stars. All images are shown with a symmetrical bi-polar stretch of ± 10 counts sec -1 pixel -1 
(F110W surface brightness of ± 2.19 mJy arcsec -2 ) to show the nature (and intensity) of the 
PSF subtraction residuals. Differential diffractive and scattering artifacts remain spatially 
fixed on the detector with the spacecraft re-orientation, but vary in intensity and phase due, 
primarily, to spacecraft "breathing" (Schneider et al. 2001). Note for example the diagonal 
diffraction spikes from the HST secondary mirror support, and a commonly occurring diffuse 
coronagraphic scattering feature (known as "the smokestack" , slightly curving to the upper 
left, dark in panels A, B & E) which appears in conjunction with a more intense "finger" 
closer to and above the coronagraphic hole (red circles), particularly bright in panels C & 
D. Panel F combines images A-D, medianed with equal weight after masking the diffraction 
spike artifacts in the individual images and rotating all to the same celestial orientation. All 
panels are 7" x 7". Green lines indicate the disk major axis. 
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Fig. 2. — A) NICMOS 1.1 two-orientation combined PSF-subtracted coronagraphic 
image of the HD 181327 debris ring (same image data as Figure IF). B) Best fit photometric 
model, as discussed in §5.4 of the text. C) Observed image (from panel a) after low-pass 
filtering with a 9° elliptical "boxcar" kernel. D) Residuals after subtracting (B) from (A). 
All panels shown at the same linear display stretch: -0.5 (black) to +10 (white) counts sec -1 
pixel -1 (-0.11 to +2.19 to mJy arcsec -2 ) and scale (as indicated). Small (r = 0'.'3) gray 
circles are centered on the location of the star. 
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Fig. 3. — ACS 0.6 /im PSF-subtracted coronagraphic image of HD 181327. A) 2"5 x 2 / .'5 
field, r = I'M (red), r = 2'.'0 (green), -0.1 to +1.5 e - s -1 pixel -1 linear display. The disk 
is reliably seen only within ± 10° of the image +Y axis (P.A. = 274?9; arcs bounded in 
orange). B) 5" x 5" field, r = I'M (red), r = 2" (green), r = 4" (yellow), to 1 e~ s" 1 pixel -1 
square root display. C) 10" x 10" field, r = 2" (red), r = 4" (green), r = 9" (yellow), [-3.3] to 
[+0.0] e~ s _1 pixel -1 log 10 display. Within ± 45° of North and 4"- 9" from the star a faint 
nebulosity can be seen, with some suggestion of curvature toward the east. 
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10 100 
Wavelength (jum) 

Fig. 4. — HD 181327 spectral energy distribution from color-corrected Spitzer/MIPS and 
IRAS photometry (see Table 2). A least-squares blackbody fit to the photometric data 
indicates a grain equilibrium temperature of 77K ± 5K. The black line is the sum of the 77K 
fit and a blackbody of T e ff = 6450K for the stellar photosphere (dotted line). 
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Fig. 5. — HD 181327 debris ring projected to a "face-on" geometry while conserving the total 
flux density. Panels (A) and (B) in this figure correspond to panel (A) and (C) in figure 2, 
but with the disk rotated to place the disk major axis (as seen on the sky) along the image 
horizontal before projection. Panel C is a linear isophotal contour map map of the image 
in panel B with 6.25% SB contour intervals from +0.1 to to +2.19 mJy arcsec -2 . Regions 
shown as hard white (SB > 2.2 mJy arcsec -2 ) are strongly affected by PSF-subtraction 
artifacts and should be ignored. The "+" symbols and small red circle mark the location of 
the star as determined from the target acquisition images. The large red circle corresponds 
to the distance at which the surface brightness peaks radially, and the blue line indicates the 
major axis of the ring as projected onto the sky. 
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Fig. 6. — Radial surface brightness profiles of the HD 181327 circumstellar ring. NICMOS 
measures along the disk minor and major axes, and in the direction of the ACS image +Y 
axis, in 0'.'5 wide strips in radial increments of 1 pixel (~ 75.8 mas). ACS measures made 
identically along the HRC image +Y axis (NICMOS and ACS images displayed at the same 
spatial scale and celestial orientation). Uncertainties in the NICMOS measures may be 
estimated from the 1 a errors shown in Figure 8. Accurately deriving the ACS measurement 
uncertainties is problematic due to the very strong azimuthal dependence of the residuals. 
Based on the profile asymmetries we estimate the uncertainties in the ACS measures in the 
1"5 - 2"2 radial region to be < ± 30% of the SB. 
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Fig. 7. — Azimuthally medianed 1.1 /iin radial surface brightness profile of the HD 181327 
ring (deprojected from an inclination angle of 31?7 from face on) and the 1 a uncertainties 
in each 75.8 pixel wide annular zone measured. The 1 a error bars indicate the dispersion 
in the NICMOS surface brightness measures around the ring at each elliptical annulus, and 
does not include the ~ 9% uncertainty in the absolute calibration of the surface brightness. 
An adjusted profile compensating for the dilution of starlight with distance from the central 
star normalized to the radius of peak brightness is also shown. 
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Fig. 8.— NICMOS and ACS deprojected SB profiles of the HD 181327 circumstellar disk 
from: Left - 2" to 4" (~ 100 to 200 AU), Right - inner detection radii to 9" (~ 450 AU). 
NICMOS profile limited to r < 4" due to field of view truncation (see figure 9) and integration 
depth. 1 a error bars are the standard deviations of the measured SB in each 1 pixel 
wide radial zone about the zonal median and do not include uncertainties in the absolute 
calibration of the disk flux densities. The ACS profile was measured after resampling the 
image to the NICMOS pixel scale (effectively smoothing the image with an 3x3 pixel 
kernel) . 



-39- 




Fig. 9.— The outer regions of the HD 181327 disk in the ACS (left) and NICMOS (right) 
PSF-subtracted images within ± 127 5 of the central star. Both images are square-root 
displays with maximum intensities (hard white) corresponding to 1.4 instrumental counts 
s _1 pixel -1 . Four stars common in both fields are noted (a-d). The debris ring revealed in 
the NICMOS image is shown at the same spatial scale and image orientation in the inset 
(linear display from -0.25 to +10 counts s _1 pixel -1 ). 
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Fig. 10. — Surface density profiles multiplied by the mean scattering cross section at 1.1 
fim assuming different vertical scale-heights (-Ho)- Left: exponential vertical profile (7 = 1). 
Right: Gaussian profile (7 = 2). 
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Fig. 11. — Midplane disk thickness in the NICMOS F110W band as a function of the distance 
from the star for different vertical scale-heights (Hq). Left: exponential vertical profile (7 = 
1). Right: Gaussian profile (7 = 2). 
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Fig. 12. — Scattered-light patterns from representative nearly-flat {(3 = 0.1) and linearly 
flaring (/3 = 1) model disks with both exponential (7 = 1) and Gaussian (7 = 2) midplane 
number density profiles for disks with vertical scale heights, H , of 1, 8, and 16 AU at r (r 
= 86.3 AU for HD 181327). Square-root stretched image displays normalized to the peak 
intensity of the ((3 = 1, 7 = 2, H = 1AU) model disk (upper right panel). 
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Fig. 13. — A representative x 2 (convergent minimization) fit of one of the MCFOST disk 
models (middle panels; /3 = 1, 7 = 2, Ho — 1 AU case as shown in Fig. 12) to the NICMOS 
observation of HD 181327 (left panels). The fitting region (in the r w l'/3 — 2'.'3 annular zone) 
is bounded by the green and yellow circles overlaid in the upper panels. The upper image 
panels are linear displays (0 to 1.5 mJy arcsec -2 , black to white, closely corresponding to the 
peak surface brightness of the HD 181327 disk. The lower image panels show the same data 
displayed over the same dynamic range but with a square-root stretch (same color-intensity 
scale as shown in Figure 12). The residuals remaining after subtracting the model image 
are shown in the linear (±1.5 mJy arcsec 2 ) display in the upper right panel. The residual 
flux densities and their 1 a deviations about their respective medians in each one pixel (75.8 
mas) wide annular zone, as a percentage of the predicted flux densities (~ 10% over the 
region of the scattered light ring), are shown in the lower right panel. 
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Fig. 14. — The minimum grain sizes and power law indices suggested, separately, from the 
observed scattered light disk asymmetry factor (§ 5.3), the [0.6] - [1.1] /im color index of the 
grains (§ 6.1), the thermal IR spectral energy distribution (§ 4), and a theoretical radiation 
force blow-out size limit (§ 8) presuming a single population of hard spherical silicate grains 
with dn(a) oc a K da. Optical constants from Weingartner & Draine (2001) are adopted for 
Mie scattering by the grains. Radiation force and grain temperatures have been calculated 
with a NextGen atmosphere model spectrum (T e ff = 6450 K, log(g) = 4.5; Hauschildt et al. 
(1999)) scaled to match the observed V-band magnitude. SED fits corresponding to reduced 
X 2 of 0.25, 1, 2.5, and 3.5 (with corresponding probabilities of 78%, 37%, 8%, and 3%) are 
shown. 




Fig. 15. — Circumstellar debris systems imaged in scattered light, (a) AU Mic, (b) (3 Pic, 
(c) HD 32297( d) HR 4796A, (e) HD181327, (f) HD 141569, (g) HD 107146, (h) a PsA. See 
Table 3 for details. 



